Dynamical Poincare Symmetry 
Realized by Field-dependent Diffeomorphisms 



00 



Oh 



00 
On 



- 1—1 

X 



R. JACKIW0 

Center for Theoretical Physics, Massachusetts Institute of Technology 
Cambridge, MA 02139-4307, USA 

A.P. POLYCHRONAKOS 

Theoretical Physics Department, Uppsala University 



D ■ S-75108 Uppsala, Sweden 

m ■ 

Faddeev Festschrift, Steklov Mathematical Institute Proceedings 



> 

' Abstract 

<N : 

We present several Galileo invariant Lagrangians, which are invariant against 
s * ' Poincare transformations defined in one higher (spatial) dimension. Thus 

these models, which arise in a variety of physical situations, provide a rep- 
resentation for a dynamical (hidden) Poincare symmetry. The action of this 
symmetry transformation on the dynamical variables is nonlinear, and in 
one case involves a peculiar field-dependent diffeomorphism. Some of our 
models are completely integrable, and we exhibit explicit solutions. 



I. INTRODUCTION 



Our colleague and friend Ludwig Faddeev has brought mathematics to physics and 
physics to mathematics. Here we mention only his work on integrable systems and 
on group theory, because these two themes also inform our essay, which is dedicated 
to him on the occasion of a significant birthday. 

We shall be concerned with Lagrangians containing fields that move in time on a 
c?-dimensional space. The models arise from diverse physical systems: continuum de- 
scription of free particle motion, isentropic fluid mechanics with irrotational velocity 
field, hydrodynamical description of quantum mechanics, free motion of membranes 
as well as higher-dimensional "<i-branes" . Our models are Galileo invariant on their 
(d, 1) space-time. However, they possess a further hidden symmetry: in terms of 
the dynamical canonical variables with which these theories are formed one can con- 
struct quantities whose algebra, determined by canonical Poisson brackets, reproduces 
a (d+1,1) Poincare algebra. Moreover, one finds that the (d + 1,1) Poincare group 
is a symmetry of the (d, 1) Galileo invariant theory, with various (d + 1, 1) Poincare 
transformations acting non-linearly on the available (d, 1) coordinates. Indeed some 



of these coordinate transformations make use of peculiar diffeomorphisms that in- 
volve the fields themselves. In summary we find that our models give a non-linear 
representation for a dynamical (hidden) Poincare group. 

Another interesting feature is that several of our models are completely integrable. 
This is seen by explicitly integrating the Euler-Lagrangian equations for some models, 
while for others - in one spatial dimension - by identifying an infinite number of con- 
servation laws, and replacing the nonlinear equations by linear ones, whose solutions 
are explicitly given. 

II. THE MODELS 

The Lagrangians that we study involve two canonically conjugate fields, 6 and p. 
When presented in first-order form, they read 

l — J d d r (op - \ P ve ■ ve - v{ P )) (2.1) 

L = J d d r (Op - \pV6 ■ V6) (2.2) 

Fields depend on time and space, {t,r); the over-dot denotes differentiation with 
respect to the temporal argument; the gradient is with respect the spatial arguments. 
Lq is the "free" Lagrangian (even though it is not quadratic), while L includes an 
interaction potential V(p), which we take to be ^-independent. 
The canonical structure implies the Poisson bracket 

{p(t,r),9(t,r')} = 5(r-r') (2.3) 

and the Hamiltonians 

H = J d d r£ 8 = \ P V9 ■ V9 + V(p) (2.4) 

H = J d d rS S = \ P V9 ■ V9 (2.5) 

(Here we are following the simplectic method for Lagrangians that are linear in the 
time derivative - this approach was advocated by Faddeev and one of the present 
authors (RJ). U) Evidently the Euler-Lagrange equations read 

p = -V • (pV9) (2.6) 

e = -\{vef + f{ P ) (2.7) 

f(p) = -j- p Jd d rV(p) (2.8) 

where the "force" term / is absent in the free case. 

We now describe several contexts that lead to this system. 
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A. Continuum Description of Free Particle Motion 

Begin with the free Lagrangian for a collection of particles, with equal masses 
scaled to unity. 

Lfree particle = | X V t^) ( m = 1 ) ( 2 ' 9 ) 

i 

Let the particle counting index i become the continuous variable r, and introduce the 
density p(t,r) and current j(t,r), with 

j = vp (2.10) 

Then ^J2 v f becomes | / d d rv 2 p = | / d d rj 2 / p. We wish to link the density and 

i 

current by a continuity equation 

p + V-j = (2.11) 

which can be enforced with help of a Lagrange multiplier 9. We thus arrive at the 
continuum Lagrangian 

£free particle ^ L = J [\] 2 / P + 9(J> + V ■ j)} (2-12) 

Since j does not participate in the canonical 1-form, / d d r6pdt = J d d r6dp, we can 
eliminate j by solving the constraint equation that follows when j is varied in L' Q [l|]. 

j = pV6 (2.13) 



Comparing (|2.10|) with (|2.13| ) shows that the velocity is given by V#; hence it is 



irrotational, with 6 playing the role of a velocity potential. 

V x v = (2.14) 

v = V0 (2.15) 
Substitution of ( gJgp in (£TJ) yields L of (gg). 



L' ^L = J d d r [9p - \pV9 ■ Vfl} (2.16) 

Supplementing Lo with a velocity-independent (^-independent) interaction term pro- 
duces @ 

L = L + J d d rV{p) (2.17) 
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B. Fluid Mechanics 



Equations ( |2.6|) , (|2.7| ) are recognized as the equations of fluid mechanics when the 
velocity field v is irrotational and the motion is isentropic. [0] The equation of motion 
(|2.6| ) is the continuity equation (|2.11|) , once ( |2.13|) is taken into account. Moreover, 



that equation (|2.7|) for the velocity potential produces the Euler equation of fluid 
mechanics 

v + v • Vv = V (pressure) (2.18) 

P 



is established by taking the gradient of Q2.7I), and using ( |2.15| ) 



v + v ■ Vv = f(p)Vp (2.19) 



The right-hand sides of ( |2.18| ) and ( [2.19] ) coincide for isentropic motion, where pres- 

~>v(p 
dp 



sure is a function only of p [Kelvin's theorem]. [In that case 9V ^ = —f{p) is the 



enthalpy w, and (p 9 gjf^ ) 2 is the sound speed u. 



C. Hydrodynamical Description of Quantum Mechanics 

For another derivation, we consider the Lagrangian for the Schrodinger equation. 

^Schrodingcr = / d d T {l^V ~ f (W)* " (W) ~ V{^)} (2-20) 

The first two terms correspond to the free, linear quantum mechanical equation, while 



V allows for possible non-linearity. Upon setting in fl2.20|) 



ip = p\e iB (2.21) 



we get || 

^schrodinger -> L = J^r {e P - \pV6 ■ V9 - V{p)} (2.22) 

where 

V{p) = V{p) + \ i ^- (2.23) 

[A similar result is obtained when a scalar field theory in (d + 1,1) dimensions is 
reduced dimensionally in the light-cone variable ~ x*- d+1 ') [0.] 
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D. d-branes 



Following Hoppe, || we consider a Nambu-type Lagrangian for a closed extended 
(i-dimensional object moving in (d + 1, 1) -dimensional Minkowski space-tine. 



L = - J d d (f)VG (2.24) 

G is (— ) d times the determinant of the induced metric 

Ox^ 1 Ox 

Ga " = Wdf? (2 ' 25) 

where a^(0) are the coordinates of the "<i-brane" in space-time (// = 0, 1, . . . , d + 1), 
depending on parameters Q (a = 0, 1, . . . , d). [0° is the evolution parameter; = 
1, . . . , d) parameterize the object at fixed time.] We introduce light-cone variables 

r =_L( x (o) + x (<i+i)) (2.26) 

6 = -L(x^ - x^) (2.27) 

and the transverse coordinates 

x i (i = l,...,d) (2.28) 

We choose the parameterization r = 0°. It follows that 

q = ( G 00 G 0s \ = ( 2d T 9 - (d T x) 2 d T 9 - d T x ■ d s x \ ^ 2g , 

al3 y G r0 —g rs J \ d r d — d r x ■ d T x —d r x ■ d s x J 

G = gT (2.30a) 

T = 2d T 6 - (d T x) 2 + g rs u r u s (2.30b) 

g = det g rs (2.30c) 

u r = d T x ■ d r x — d r 9 (2.30d) 

Here d T = ^ = d r = r = 1, . . . , d; g rs is the inverse of g rs = d r x ■ d s x and will 
be used to move the (r, s) indices. Note that the dimensions of the x space (indexed 
by i) and of the parameter space (indexed by r) are both d. The Euler-Lagrange 
equations are conveniently presented in canonical form 

d T x = -p/n + u r d r x (2.31) 
dr0=^(p 2 + g)+u r d r e (2.32) 

d T p = -d r (^gg rs d s x) + d r (u r p) (2.33) 
<9 T n = d r (Uu r ) (2.34) 



with II and p satisfying the constraint 



p ■ d r x + Ud r 9 = (2.35) 

Here u r is as given by ( |2.30d| ) [this is a consequence of ( |2.31| ), (|2.35|) 1, and can be set 
to zero by appropriate choice of the parameterization. Thereupon it follows that d T Tl 
vanishes, so we choose II to be — 1. The equations then reduce to 

d T x = p (2.36) 
d T e=\{p 2 + g) (2.37) 
d T p = d r (gg rs d s x) (2.38) 



and the constraint reads 



d r 6 = p-d r x (2.39) 



The constraint is now solved by a transformation introduced by Bordemann and 
Hoppe. Rather than viewing the (dependent) variables p and 9 as functions of the 
(independent) parameters r and 0, a change of variables is performed by inverting 
x(t,4>) and expressing <p in terms of r and x, (renamed r), so that p and 9 be- 
come functions of r (renamed t) and r. It then follows from the chain rule that the 
constraint ( |2.39| ) becomes 

dr l 89 _ i dr % 

and is solved by 

p = V9 (2.41) 

where the gradient is with respect to r. Moreover, according to the implicit function 
theorem (see sidebar) 

d T = d t + V9 ■ V (2.42) 

so that (|2.36|) reproduces Q2.41|) , since d t x l = 0; while (|2.37| ) and (|2.38[) read, respec- 
tively 

9 + \V9-V9 = \g (2.43) 



p + V9 ■ Vp = V(0 + \V9 ■ V9) = d r (gg rs d s r) 



(2.44) 



The over dot signifies differentiation with respect to t. 

It remains to show that f^43"D and ( $M ) are consistent with each other. This is 
achieved by recognizing that 
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nn rs — / ,ri 2 ...i d sj 2 -jd n „. . (n A^ a \ 



TI) J Ox^ Ox^ 

OUt Qij = -QTiQTJ, SO 



„„»f) ^ - 1 r ri2...i, dxk2 ^ r s.h...M dxi ^ dxkd 

yy — . . t . . . t 



(d- 1)! ••• dft* d(p s d<j)K 



r, 2 .., d ^_ . . . ^_^...k d ^ ^ (2 45b) 



(d-1)! ^ 



i 



Note that det fe- = g 2 . It now follows that 



2 c « 2 ...i d jk 2 ...k d dxk2 . . . dx k *dgi dx* 



(d-l)\ d<f>^ d<fh dxi 



1 

I 1 (9(7 2 

Jk 2 ...k d Jk 2 ...k d g2 J 



{d-l)\ y 9xi 

1 dg 

2 dx % 



(2.45c) 



Thus ( p.44j ) implies ( |2.43| ), once an r-independent constant of integration is absorbed 
in 6. 

It is seen that the equations for the "(i-brane" collapse into ( |2.43|) . But we still 
need an equation for g, which can be obtained by differentiating g with respect to r 
and using ( |2.42|) . 



gg rs -^g rs = d T g = g + V6-Vg (2.46) 



But 



d 

9 rs 7r9rs = 2g rs d r x l d T d s x l 
or 

= 2g rs d r x i d s p i 

Since g rs = f£g£, it follows that the right side of ( ^46|) is 2g^d s p i = 2g£ T p i = 
2gV 2 9. Therefore if we define g = 2X/p 2 , (|2.46|) is equivalent to 

p + V • (pVG) = (2.47) 

We conclude therefore that the motion of a "d-brane" , moving in time on d + 1 
spatial dimensions is governed by equations derivable from our Lagrangian (|2.1|) , with 

V(p) = X/p. (2.48) 
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Sidebar: Derivation of (|2.42|) Consider a function / that depends on r and 
<fi: f = /(r, <j>). The r derivative differentiates / with respect to the first argument. 
Next express as a function of r and x, renamed r. It is clear that 



d T f(r, cp) = ^/(r, 0(r, r)) - ^/(r, 0(r, r))^|^ 



(2.49) 



Next observe that 



When is held constant, we have 



d<f> d<p dr l 
dr dr 1 dr 



(2.50) 



(2.51) 



constant 



But from fl2.36|) it follows that the r derivative of r is p, which according to (|2.41 ) is 
V# . Thus ( p.51|) becomes 



30 _dcj)_d6_ 

dr dr 1 dr 1 



(2.52) 



Substitution in (|2.49|) and use of the chain rule shows that (|2.42|) is true, once r is 
renamed t. 



d „ d6 d(j) s df 



r=t dt dr* dr 1 d(p s 



(2.53) 



III. SYMMETRIES OF THE MODEL 



A. Galileo Symmetries 



From the derivations, it should be obvious that Lq ( |2.2j ), and also L Q2.1| ) (with 
obvious restriction on V) possess the Galileo symmetry. Let us record the generators 
of the infinitesimal transformations as integrals of the appropriate densities; also we 
specify the action of the finite transformations (parameterized by to) on the fields: 
p — > p w ,Q — > 6 W , by presenting formulas for p u (t,r) and ^(t, r) in terms of p(t,r) 
and 6(t,r). One verifies that the generators are time independent according to the 
equations of motion (p.6|) - (p.8[) , and this further implies that the transformed fields 
p w and U also solve ( p.6|) -( p78|) , when p and 6 are solutions. 
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Time, space translation 
— Energy 

H = Jd d rS, £ = \pV9-V9 + V(p) = \j 2 /p + V(p) (3.1) 



Momentum 



P = Jd d rV, V = pV6=i (3.2) 



Space rotation 

— Angular momentum 

J*J = / d d r [r l V 3 - r 3 V l ) (3.3) 



With these space-time transformations p u , 9^ are obtained from p, 9 by respectively 
translating the time, space arguments and by rotating the spatial argument. f| 

• Galileo boost 



Boost generator 



The boosted fields are 



B = tP 



Jd d rrp (3.4a) 



p u (t,r)=p(t,r-ort) (3.4b) 
M*> r ) = Ofa r-u>t)+u>-r- ujH/2 (3.4c) 
The inhomogeneous terms in 9 U are recognized as the well-known Galileo 1-cocycle, 



compare ( |2.21| ). Also they ensure that the transformation law for v = V# 

v(t, r) -> v w (t, r) = v(t, Y-uot)+u (3.4d) 

is appropriate for a co-moving velocity. Furthermore, knowledge about the Galileo 
2-cocycle leads us to examine the P, B bracket, and its extension exposes another 
conserved generator, arising from an invariance against translating 9 by a constant; 



this just reflects the phase arbitrariness in (|2.21| ). 
• Phase symmetry 

— Charge 

N = J d d rp (3.5a) 

Pw= P (3.5b) 

9 u = 9-uj (3.5c) 
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B. Connection with Poincare Symmetry 



It is well known that a Poincare group in (d + 1, 1) dimensions possesses the above 
extended Galileo group as a subgroup. This is seen by identifying selected light-cone 
components of the Poincare generators P^, with the Galileo generators, 

p m = ( P - ) P + ) P i) ~ (h } N, P l ) (3.6) 

= (M + ~,MT\M + \M ij ) 
M +l « B\ M ij « J*' (3.7) 

where the ± components of tensors are defined as in ( |2.26| )-( f2.27| ) 

T (±) = _L ( T (o) ±T (<m)) . (3. 8 ) 

(This fact is responsible for behavior in the "infinite momentum" frame.) |§ But the 
Lorentz generators M + ~ and M~ l have no Galilean counterparts. 

A remarkable fact, first observed in Ref. and then again in simplified form in 
Ref. |J, is that in the model ( [2.1|) with V(p) = A/p [as in the "<i-brane" application, 
eq. (|2.48|) 1 one can define quantities that can be set equal to the generators missing 



from the identification ( |3.7| ), namely M + ~ and MT\ This holds for arbitrary interac- 
tion strength A; setting it to zero allows the same construction for the free Lagrangian 
(|2.2j ). We wish to determine what role the additional generators have for the models 



(2T) and Q. [10| 



C. Additional Symmetries 

We observe that the free action I = J dtL , as well as the interacting one 

h = J dt d d r (Op - \pV6 ■ V6 - X/p) (3.9) 

are invariant against time rescaling t — > e^t, which is generated by 

D = tH-Jd d rp6. (3.10a) 

Fields transform according to 

p(t, r) -> Pu (t, r) = e~"p{e% r) (3.10b) 

6(t, r) -> 9 u (t, r) = e^e^t, r) . (3.10c) 

The dilation generator D is identified with M + ~. It is straight forward to verify from 
(|2.6|)-(|2~8|) that D is indeed time independent. 
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More intricate is a further, obscure symmetry whose generator can be identified 
with M~\ Consider the field-dependent coordinate transformations, implicitly defined 
by 



t -> T(t, r) = t + uj-r + \u 2 6(T, R 
r -> R(t,r) = r + u>6{T,R) 



(3.11) 



with Jacobian |J|. 



J = det 



/ dT dT \ 
~dt ~dr~i 
dR 1 dR i 



l-u- V0(T, R) - W6(T, R 



(3.12) 



The transformation parameter uj has dimensions of inverse velocity. When fields are 
taken to transform according to 

1 



p(t,r) ^ Pu! (t,r) = p(T,R) — 



(3.13a) 

e(t,r)^e w (t,r) = e(T,R) (3.13b) 
one verifies that I\ and Iq are invariant. This is readily seen for the interaction term 



A / dtd d r 



P(t,r) 



A / dt d d r 



P(T, R) 



A / dTd d R 



P(T, R) 



(3.14) 



To establish invariance of In, it is useful to write it first as 



J dtd d r p(9 + \V9 ■ V0) 



(3.15) 



The desired result follows once it is realized that the quantity in curly brackets equals 
J 2 {9(T,R) + |(V0(T,i?)) 2 }. The transformations ( |3TTB| ) are generated by 



G 



J d d r (r£ - \pV6 2 
J d d r {rS - 6V) 



(3.16) 



which is time independent according to ( |2.6j )-( |2lj| ), and whose algebra with the other 
generators show that G l can be identified with M~\ 

While we have no insight about the geometric aspects to this peculiar symmetry, 
the following remarks may help achieve some transparency. 

Observe that the Galileo generators can be expressed in terms of p and j or 
p and v = j jp — V0. Consequently, they are also defined for velocity fields with 
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vorticity, (Vx v 7^ 0), and provide well-known constants of motion for the (isentropic) 
Euler equations [3]. However, the velocity potential 9 is needed to form D and G, 
which therefore have a role only in vortex-free motion (with a specific potential or no 
potential). 

From the identification with Poincare generators, we see from ( |3.4| ), ( |3.6| ), ( |3.1(J| ) 



and ( |3.16|) that p is the P + density and that 9 plays the role of x . The latter 



identification is further suggested by the fact that for all the transformations that 
are identified with Lorentz transformations, viz. (|3.3|), (|3.4|), ( |3.10| ), ( |3.11| ), ( |3.13| ) it 
is true that 

2T9(T, R)-R 2 = 2t9 u {t, r) - r 2 

where T and R are the appropriately transformed coordinates. The naturalness of 
this expression is appreciated when it is recognized that 

In the "(i-brane" development we saw that 9 is indeed x~, see ( |2.27|) . Moreover 
there - oc yfg is the Jacobian of the inversion transformation (r, <fi) — » (t, r(t, <fi)), so it 

is quite natural that under the further transformation (£, r) — > (T(t, r), R(t, r)\ 1/p 
acquires the Jacobian of that transformation. Presumably transformations (|3.11|) , 



( 3.1 3| ) reflect a residual invariance or gauge- fixed "<i-brane" theory, but the reason for 



the specific form ( |3.11|) of the transformation is not apparent. 



In applications to fluid mechanics our transformation generates nontrivial solu- 
tions of Euler's equations, as we now explain. 

IV. EXPLICIT SOLUTIONS AND THEIR TRANSFORMS 



In order to gain insight into the peculiar diffeomorphism transformations ( 3.1 1| ) 



(|3.13| ), we consider its effect on some explicit solutions to Eqs. ( |2.6|) -( [2~"8D , in the free 
{V = 0) and interacting (V = X/ p) cases. 

A. No Interaction, V = 

1. Particular Example 

With V = 0, eq. ( ^.7| ) decouples from (2.6), and is solved by 

2 

9(t,r) = fj (4.1a) 

which, apart from selecting an origin in time and space and presenting a rotation and 
boost invariant profile, is also invariant against time rescaling ( |3.10| ) and the uncon- 
ventional diffeomorphism ( |3.11| )- (|3.13j ). The fluid moves with a velocity unaffected 



by boosts (|3T4H|) 
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r 

v = — 

t 



(4.1b) 



The density is not determined, since the solution of the continuity equation ( |2.6|) in d 
spatial dimensions involves an arbitrary function of t/r, and of the angles specifying 
r that are denoted by f = r jr . 



P(t,r) 



f(t/r,r) 



(4.1c) 



With 9 as in ( |4.1a| ), the coordinate transformations ( |3.11| ) takes explicit form 



T(t,r) 
R(t,r) 



J 



tr l 



rr„ 



r 



u r 
r H 

2 t 



(4.2) 



so that, as stated, the transformed 8 w (t, r) = 9(T, R) = R 2 /2T coincides with 9(t, r), 
while the transformed density retains the form ( |4.1oj ) but with a different function of 

t/r, 



f(tr u /r, 



(4.3) 



This coincides with p(t,r) for the special choice f(t/r,r) oc (t/r) + in ( 4.1c ), which 



provides a density profile that is invariant under the diffeomorphism (|3.11| ' 

One may construct other "free" solutions, for which 9 remains unchanged under 



(HI). 



( |3.11| ), ( |3.12| ), ( p,13b|) , while p involves arbitrary functions. Also there are free so- 
lutions that respond nontrivially to the transformations. We do not pursue specific 
solutions any further, because it is possible to give the general solution to the free 
problem, which we now present. 



2. General Solution 

Since the pair of equations ( |2.6| ), (|2.7|) is first-order in time, to give a general 
solution we need initial data at initial time t Q ; owing to time-translation invariance, 
t D can be taken to be 0, without loss of generality. 

p(0,r) = p o (r) (4.4) 
0(O,r) = o (r) (4.5) 

The general solution to the free version of fl2.6Q , Q2.71) , that is, at / = 0, is given in 
terms of a vector- valued function of t and r, q(t,r), which satisfies 

13 



q + tV9 (q) = r 



(4.6) 



The solution to (|2.6|), (|2.7|) then reads 



p(t,r) = p (q) 



det 



9(t,r) = 9 (q) + ~(y9(qf 



(4.7) 
(4.8) 



Note that the velocity, which is the gradient of 9, is just the initial velocity, evaluated 
on q. 



v(t,r) = V r 9(t,r) = V q 9 Q (q) = v G (q) 



Consequently, ( |4.6| ) may also be presented as 



q + tv (q) = r 



(4.9) 



(4.10) 



and one verifies that the free Euler equation, without the irrotational condition on v, 
Vxv/0, 



v + v Vv = 



(4.11) 



is solved by v (q), where v (r) is the initial velocity and q satisfies ( |4.1(J| ). This of 
course is nothing but the description of freely moving dust particles. Also one verifies 
that the form ( (4.7[ ) for p solves the continuity equation, even when V x v ^ 0. 
Note that the free motion of dust particles leads to the conserved quantities 



J f = J d d rp(t,r)f(v«{t,r), 
j f = 



v id (t,rj) 



(4.12) 



where / is an arbitrary function of the velocity components. This is an expression in 
the continuum formalism of the fact that velocities of the free theory are constant. 
At d — 1, there are additional conserved quantities 



dx f(v(t, x)) 



(4.13) 



for arbitrary functions /. 



14 



B. With Interaction, V = X/p 



1. Particular Example 

With interaction, the equations of motion remain coupled. 

p = -V ■ (pV9) (4.14) 

e = -\ve- ve + x/ P 2 (4.15) 

Remarkably a solution exists that is similar to (|4.1|) : demanding rotational and time- 
rescaling invariance, [4.14| ) leads to a second order equation in r for 9, once p has been 



eliminated with the help of ( 4.15|) . A particular solution of that equation then gives 



In \ \ + \ 

p(t,r) = y-j(d-l)^ (4.17) 
The solution is not diffeomorphism invariant; see below. The velocity flow is 

(4.18) 



(d-l)t 



while the current reads 



/2A. , 

where the sign is determined by the sign of t. Evidently here d must be greater 
than 1 and A must be positive. (The latter requirement is natural, since in the fluid 
mechanical application the sound speed u for our model is = y/2X/p.) 



2. General Solutions at d = 1 
It turns out that this model, along with a much more general class of models 



with local interactions, is completely integrable in one dimension. [H] The integrals 
of motion can be expressed in the form 

j(±) 



J dxpf ± 



if> = (4.20) 
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for arbitrary functions /. We obtain two 'chiral' sectors involving the velocity of the 
fluid plus (minus) the local velocity of sound. 

The general solution can be obtained in this case through linearization. It is 
known the non-linear Euler and conservation equations can be mapped into a linear 
system. [^] This is achieved by performing a double Legendre transform on 9, thereby 
replacing the independent variables t and x, with v — |^ and with the enthalpy 



w 



_de _ if §e 

dt 2\dx 



dx 
2 



x „, , d9(t,x) dO(t,x 
yfa w ) = 9(t, x) - t — - x 



dt dx 
9(t,x)+t(w + \v 2 )-xv (4.21) 



It is then true that 



dw dw dv 

An equation for ^ follows from the continuity equation for p, when it is remembered 
that for isentropic motion p can be taken to be a function of w. In our model 
w = —X/p 2 . One finds 

u2 d 2 V d 2 ^ + d^> _ Q ^ 4 23 ^ 

dw 2 dv 2 dw 

where the sound speed u = (p^-) 1 ^ 2 is regarded as a function of w. For our problem 
u 2 = 2X/p 2 = —2w. Note that the coupling strength A does not occur in (f4.23|) . It 
reenters the formalism through w = —X/p 2 , w ^ for A ^0. 

The general solution of ( |4.23| ) can be expressed in terms of two general functions, 
F and G, of one variable. It is straightforward to verify that the following expression 
solves Q4.23Q , with u 2 = -2w. 



w) = F(v + y/-2w) - V-2w F'{y + y/-2w) 

+ G(v - v 73 ^) + V^2w~ G'{v - V^2w~) (4.24) 

Explicit solutions can be obtained from this expression for specific choices of F, G. 

We record the general time-rescaling invariant solutions. The time-rescaling 
Ansatz for 9 and p [9 oc 1/t, p oc t] lets us use ( |4.15| ) to express p in terms of 9 
and its spatial derivative. Then ( [4.14| ) yields a second-order equation in x for 9, 
which is solved by expressions that involve two arbitrary integration constants. For 
A > we find 

1 I 

9(t,x) = —— sinh 2 fcc, --— cosh 2 A;x (4.25a), (4.25b) 



X 2A A: |t | V2Afc|t| 
sinh 2 kx ' cosh 2 kx 
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= ■ ■ ___ 1 2,'.l (4.26a), (4.26b) 



while for A < 



6 ^ x ) = 7^7^ 2 kx, -^cos 2 kx (4.27a), (4.27b) 



2\X\k\t\ J2\X\k\t 



= VfcT' ^kT- (4.28a), (4.28b) 

Here k is an arbitrary, positive integration constant; the second constant gives an 
origin to x, and has been suppressed in the above by setting it to zero. Note that the 
current of the second solution at A > has the kink shape. 



j = =Fv2Atanh fcx (4.29) 

This profile puts one in mind of soliton phenomena, and hints at the integrability of 
Eqs. (|4.14j ), (|4.15|) in one dimension. 



In terms of the linearized formalism of ( |4.21| )- (|4.24| ), the A > solutions (4.25a) 
and (4.26a) correspond to 

F{z) = -G{z) = ^{l-\nz) (4.30) 
while (4.25b) and (4.26b) have 

F'(z)=G(-z) = ^(l-\nz) (4.31) 
The A < solution in (4.27) and (4.28) are the analytic continuation of the above. 



3. Transforming the d = 2 Solution 

We now exhibit the form of the solutions when the field-dependent, coordinate 
transformation (|3.11| )- (|3.13|) are performed on ( |4.16|) -( |4TT9"1) . For simplicity we discuss 



only the d = 2 (membrane) case, and take t > 0. The new coordinates are determined 
by the old ones by ( PH) and (|4l6|) . 



T = \t + \u ■ r ± \yj{t + 2u> • r) 2 - 2u 2 r 2 (4.32a) 
1 1 lj-R uo 2 R 2 , Ann ^ 



R = r + — \-t - 2u ■ r ± yj(t + 2u ■ r 2 - 2u 2 r 2 } (4.32b) 



After 6 and p are transformed according to the rules (|3.13j ), it is noticed that ex- 
pressions are simplified by performing the Galileo boost r — > r — ujt/u 2 , according 
to the rules (|3.4j) . (This precludes taking the limit u) — > 0.) Also, time is rescaled 



according to (|3.10|) , with t — > \/2t. Finally, we define lj/cu 2 = c, which has dimension 
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of velocity, and then the transformed profiles provide two solutions, depending on the 
sign of the square root 



9 c (t, r) = ±^2{c ■ r) 2 - c 2 r 2 - cH 2 

2(c • r) 2 - c 2 r 2 - cH 2 



, s V2A 



The velocity is 



v c (t,r) 



r 2 =f 2tJ2(c • r) 2 - c 2 r 2 - c*t 2 



2c(c • r) — rc 2 
2(c ■ r) 2 - c 2 r 2 - cH 2 



1/2 



(4.33a) 
(4.33b) 



(4.33c) 



and the current reads 

j c (t,r) = ±V2A- 



26(6 • r) — r 



[r 2 =F 2t J2(c ■ r) 2 - c 2 r 2 - cH 2 } 1 / 2 



(4.33d) 



Note that c may be replaced by ic and p c by — p c , to obtain another solution. 

In the figures we exhibit the profiles of the interacting solutions. We plot the 
original and transformed densities, and the transformed currents j c = p c V8 c , in terms 
of the variables r/t (t > 0). Without loss of generality, c is taken along the x-axis, and 
its magnitude is incorporated in the dimensionless ratio r/ct. The original density 
possesses a singularity at the origin; in the transformed solutions the singularity is 
present only with the upper (negative) sign in the bracketed expression of ( [4.33b| ). 
where its denominator vanishes at r 2 = (c ■ r) 2 = 2c 2 t 2 . The transformed currents 
exhibit a similar singularity. In the physical region the argument of the square root 



must be positive, 2(c • r) — r - 
validity for some of the profiles. 



crt > 0. This requirement creates an envelope of 
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FIG. 4. The transformed current j c (t,f)/V2A, with the upper signs. The envelope 
defining the physical region is at x 2 — y 2 = c 2 t 2 . 
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FIG. 5. The transformed current j c (i, r)/v2A, with the lower signs. The envelope defin- 
ing the physical region is at x 2 — y 2 = c 2 t 2 . 
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